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Abstract 


V 


\  Accurate  radar  tracking  of  targets  flying  at  low  altitudes  above  a  smooth 
surface  is  difficult  because  of  the  surface  reflection.  We  propose  a  solution  based  on 
deterministic  physical  modelling  of  the  specular  multipath  and  the  maximum 
likelihood  method.  This  paper  describes  the  techniques  and  the  results  of  performance 
studies.  We  derive  the  Cramer— Rao  bound  and  show  the  benefit  of  using  the  refined 
propagation  model.  Monte-Carlo  simulations  are  employed  to  compare  the 
performance  vdth  the  Cramer-Rao  bound  and  demonstrate  threshold  effects  on  target 
height  estimation.  Finally  we  present  the  results  obtained  with  two  X— band  radar 
experimental  systems.  _ 

Ri:suM^ 


Le  pistage  precis  d’une  dble  radar  volant  au-dessus  d’un  plan  d’eau  est  rendu 
difficile  par  la  presence  des  reflections  sur  la  surface.  Nous  proposons  dans  ce  cas 
d’indure  un  modde  speculaire  de  la  propagation  multivoie.  Ce  modde  est  utilise  avec 
le  Maximum  de  Vraisemblance.  Dans  ce  rapport  nous  presentons  la  technique, 
derivons  la  limite  de  Cramer-Rao  et  pr^entons  des  r^ultats  de  simulation  indiquant 
que  I’emploi  d’un  modde  speculaire  ameliore  grandement  la  precision.  Nous 
pr4sentons  finalement  des  r^sultats  exp^rimentaux  obtenus  par  deux  systemes  radar 
operant  dans  la  bande  X.  Nous  demontrons  que  I’utilisation  d’un  modde  speculaire 
donne  un  bon  pistage  lit  ou  I’approche  conventionelle  ne  fonctionne  pas. 
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EXECUTIVE  SUMMARY 


This  report  presents  a  new  method  for  low— angle  radar  tracking  in  the  presence 
of  specular  multipath.  We  define  new  estimators  based  on  the  maximum  likelihood 
criterion  for  estimating  the  height  of  targets  flying  at  low  altitudes  over  the  sea.  We 
assume  that  interfering  multipath  signals  are  present  and  show  that  identifying  and 
modelling  the  nature  of  the  multipath  interference  can  have  a  beneficial  effect  on  the 
estimation  of  desired  parameters. 

Previous  experimental  works  indicated  that  at  low  altitude  the  interference 
pattern  can  be  accurately  predicted  under  normal  propagation  conditions.  We  use  that 
information  in  an  optimal  way  in  the  context  of  the  maximum  likelihood  (ML) 
estimation  theory  and  we  name  this  algorithm  the  Refined  Maximum  Likelihood 
(RML)  technique. 

We  use  real  multipath  data  to  test  the  RML  technique.  In  recent  years,  DREO 
(Defence  Research  Establishment  Ottawa,  Canada)  has  been  undertaken  many  trials  to 
collect  real  multipath  data  on  the  Ottawa  River.  As  part  of  the  ELAT  (Experimental 
Low— Angle  Tracking  radar),  an  8-element  sampled  aperture  antenna  was  used  to 
collect  the  data  at  two  frequencies:  8.6  and  9.6  GHz.  The  target  was  a  high-power 
beacon  source.  Also,  in  late  October  1987,  a  series  of  experiments  was  performed  on  the 
west  coast  of  the  Bruce  Peninsula  overlooking  Lake  Huron  (Canada).  This  time,  the 
MARS  (Multi-parameter  Adaptive  Radar  System)  developed  at  the  Communications 
Research  Laboratory  (CRL),  McMaster  University  (Canada)  was  used.  MARS  has  a 
32-element  sampled  aperture  antenna  which  was  used  to  collect  data  at  multiple 
frequencies  within  the  8-12  GHz  band.  Here  again,  the  target  was  simulated  by  a 
beacon  source.  Data  collected  by  these  two  systems  form  our  experimental  data  base  to 
test  our  technique. 

The  results  of  simulations  and  experiments  indicate  substantial  improvements 
with  the  RML  algorithm  compared  with  conventional  techniques  such  as  Fourier 
beamforming.  The  results  show  some  robustness  of  the  RML  technique  in  the  presence 
of  diffuse  multipath  with  experimental  results  obtained  in  sea  state  3  conditions. 
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A  RefiBed  Maximnm  LUcdihood  Method  for  Tracking  Low— Altitude 
Targets  ovk  the  Sea:  Results  of  Simulation  and  Experiments. 


1.0  INTRODUCTION 


This  paper  presents  a  new  method  for  low— angle  tracking  in  the  presence  of  specular 
multipath.  We  define  new  estimators  based  on  the  maximum  likelihood  criterion  for 
estimating  the  height  of  targets  flying  at  low  altitudes  over  the  sea.  The  problem 
investigated  is  direction  finding  for  a  single  target  in  the  presence  of  specular  reflection 
from  a  relatively  smooth  flat  smface  such  as  the  sea.  In  recent  years  a  considerable  amount 
of  research  has  been  done  in  the  area  of  signal  processing,  using  superresolution  algorithms 
to  separate  a  target  from  its  image  and  to  provide  an  accurate  estimate  of  the  target’s 
elevation.  Our  approach  is  essentially  different;  we  assume  that  interfering  multipath 
signals  are  present  and  show  that  identifying  and  modelling  the  nature  of  the  multipath 
interference  can  have  a  beneficial  effect  on  the  estimation  of  desired  parameters. 

Barton  [1]  summarizes  the  low-angle  tracking  work  up  to  1974.  His  paper  presents  a 
number  of  methods  to  combat  miiltipath  errors.  Most  of  these  methods  use  monopulse  and 
they  all  fail  to  resolve  targets  separated  by  less  than  0.25  beamwidth.  More  recently 
Haykin  and  Kesler  [2,3]  propose  an  adaptive  canceller  and  Haykin  [4]  investigates  the  use 
of  high  resolution  techniques  such  as  autoregression,  maximum  entropy  and  linear 
prediction.  It  seems  [5]  that  a  limit  of  0.25  beamwidth  applies  to  most  high  resolution 
techniques  [6],  including  the  well  known  Multiple  Signal  Classification  (MUSIC) 
method  [7].  MUSIC  and  most  of  these  techniques  fail  to  resolve  correlated  signals  such  as 
those  encountered  with  multipath;  here  the  direct  ray  is  fully  correlated  with  the  reflected 
ray  and  the  required  resolution  can  be  as  high  as  than  0.05  beamwidth. 

There  are  a  number  of  more  or  less  efficient  techniques  to  decorrelate  the  signals. 
Filial  [8]  discusses  techniques  using  spatial  smoothing.  With  this  technique,  we  form 
sub— arrays  and  average  the  corresponding  covariance  matrices.  Filial  shows  that  to  detect 
the  arrival  angles  for  two  closely  spaced,  equi— power,  correlated  signals,  requires  increasing 
the  number  of  sample  vectors  (snapshots)  by  a  factor  of  approximately  [(1/Kwd)*-1]  in 
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comparison  with  the  number  required  for  uncorrelated  signals.  Parameter  K  represents  the 
number  of  sensors  and  <Jd  =  M^cos  0i  —  cos  0^j2  is  related  to  angular  separation  between 
signals.  For  low— angle  tracking,  we  would  require  a  considerably  increased  number  of 
snapshots.  The  factor  is  larger  than  10000  since  the  angular  separation  between  the  target 
direction  9\  and  the  image  direction  62  may  be  less  than  one  degree.  Other  techniques  [9,10] 
use  a  multidimensional  search  vrith  multiple  parameters,  imposing  a  high  computational 
burden  that  may  be  impractical  for  low— angle  tracking. 

We  propose  a  different  technique  which  incorporates  a  highly  descriptive  model  of 
the  physics  of  the  low— angle  problem.  Previous  experimental  works  [11—13]  indicated  that 
at  low  altitude  the  interference  pattern  can  be  accurately  predicted  under  normal 
propagation  conditions.  Litva  [14—16]  first  used  a  similar  model  in  an  algorithm  called 
CHA  (Correlation  Height  Analysis)  and  obtained  resolution  up  to  1/28  of  a  beamwidth  [16] 
under  normal  propagation  conditions  and  smooth  sea  surface.  We  use  the  idea  of  a  detailed 
propagation  model  but  in  the  context  of  the  maximum  likelihood  (ML)  estimation  theory 
and  we  name  this  algorithm  the  Refined  Maximum  Likelihood  (RML)  technique  [17—19]. 

In  this  paper,  we  present  the  first  development  phase  of  the  RML  technique  where 
the  propagation  model  assumes  a  medium  that  is  linear,  homogeneous,  isotropic  and 
frequency  invariant  (narrow— band).  The  effects  of  anomalous  propagation,  ducting,  sea 
roughness  and  sea  swell  are  under  investigation  and  results  will  be  reported  as  a  second 
paper.  The  RML  technique  is  designed  to  work  in  situation  where  other  high  resolution 
methods  fail  to  resolve  the  target  from  its  image.  If  the  angular  separation  between  the 
target  and  its  image  is  greater  than  one  beamwidth,  a  monopulse  or  Fourier  beamforming 
method  can  be  used.  This  is  simpler,  more  robust  and  computationally  more  efBcient  than 
the  RML  technique. 

Maximum  likelihood  estimation  has  been  applied  to  angle-of-arrival  estimation 
using  an  angle— of— arrival  (AOA)  model.  In  [20],  a  comparative  study  of  a  number  of 
apparently  different  angle  estimation  techniques  [21—28]  showed  that  they  were,  in  fact, 
ML  techniques  using  the  AOA  model. 

The  AOA  model  has  an  unknown  angle  and  an  unknown  complex  amplitude  for 
each  signal;  specular  multipath  is  represented  by  two  such  signals  in  the  AOA  model 
yielding  six  parameters  which  must  be  implicitly  estimated.  In  contrast,  we  use  a  priori 
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infonnation  to  reduce  the  number  of  unknowns  in  our  propagation  model;  specular 
multipath  is  modelled  with  only  three  unknowns:  an  amplitude,  a  phase  and  the  target 
height.  To  obtain  the  information  required  by  the  model,  the  radar  is  operated  in  an 
acquisition  mode  to  determine  target  range  and  doppler  velocity.  Prior  knowledge  about 
the  geometry  and  the  specular  reflection  coefficient  is  assumed.  This  is  not  a  very 
restrictive  assumption  for  low— angle  tracking  because  the  reflection  coefficient  is  almost  —1 
under  smooth  sea  conditions.  Sea  roughness  and  divergence  due  to  surface  curvature  are 
taken  into  account  by  decreasing  the  amplitude  of  the  specular  reflected  ray.  Using  this 
model  we  have  accurately  tracked  beacon  targets  over  rough  sea  conditions  (sea  state  3). 

We  use  real  multipath  data  to  test  the  RML  technique.  In  recent  years,  DREO 
(Defence  Research  Establishment  Ottawa,  Canada)  has  been  undertaken  many  trials  to 
collect  real  multipath  data  on  the  Ottawa  River.  As  part  of  the  ELAT  (Experimental 
Low— Angle  Tracking  radar),  an  8-element  sampled  aperture  antenna  was  used  to  collect 
the  data.  The  target  was  a  high-power  beacon  source.  Also,  in  late  October  1987,  a  series 
of  experiment  were  performed  on  the  west  coast  of  the  Bruce  Peninsula  overlooking  Lake 
Huron  (Canada).  This  time,  the  MARS  (Multi— parameter  Adaptive  Radar  System) 
developed  at  the  Communications  Research  Laboratory  (CRL),  McMaster  University 
(Canada)  was  used.  MARS  has  a  32-element  sampled  aperture  antenna  which  was  used  to 
collect  data  at  multiple  frequencies.  Here  again,  the  target  was  simulated  by  a  CW  beacon. 
Data  collected  by  these  two  systems  form  our  experimental  data  base  to  test  our  technique. 

The  paper  is  organized  as  follows.  Section  2  presents  the  propagation  model.  In 
Section  3,  we  derive  the  RML  height  estimators.  Section  4  contains  the  derivations  of  the 
Cramer— Rao  bounds.  Section  5  gives  the  results  of  simulations  and  Section  6  terminates 
with  the  experimental  results.  We  use  the  following  notation:  matrices  are  represented  by 
bold  upper-case  letters,  vectors  by  bold  lower-case  letter,  scalars  by  both  upper  and  lower 

italic  letters.  The  superscripts,  **,  ||  ||  denote  estimate,  conjugation,  transposition, 

conjugation— transposition,  and  vector  norm  respectively.  Also  signifies  a  M>«N 
complex  matrix. 


2  The  Propagation  model 


The  propagation  model  used  in  this  paper  assumes  a  medium  that  is  linear, 
homogeneous,  isotropic  and  frequency  invariant  (narrow— band).  The  model  considers 
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bipath  propagation  with  only  a  single  reflected  ray  emanating  from  a  virtual  target  image. 
We  assume  that  the  targets  are  located  at  such  a  distance  from  the  receiving  antenna  that 
the  impinging  waves  can  be  considered  as  being  planar.  We  consider  a  two-way 
transmission  model  (radar  model)  and  a  one-way  transmission  model  (beacon  model).  The 
beacon  model  is  important  for  experimental  investigations.  It  is  much  easier  and  less 
expensive  to  set  up  a  beacon  system  and  a  receiving  array  with  narrow— band  receivers  than 
to  implement  a  radar  system  vrith  its  much  more  complicated  transmission  and  reception 
system.  The  experimental  results  presented  in  this  paper  are  obtained  using  a  beacon  as  a 
target;  experimental  verification  using  a  radar  transmitter  is  ongoing. 

The  geometry  for  the  beacon  model  is  illustrated  in  Fig.l  assuming  an  equivalent 
flat-earth  model  with  parameters  compensated  for  the  effects  of  the  earth’s  curvature.  The 
noise— free  observation  model  for  a  signal  received  at  the  element  of  an  array  of  K 
sensors  is 

=  Qi  exp(-jiR)  exp{-j({z^  +  hl)/2R}  ^ 

[exp(j(h\,Z[,/R)  +  A  exp(-j(hi,Zk/R)],  (2.1) 

where  ^  =  2ir/X,  X  is  the  wavelength,  Zk  is  the  height  of  the  element,  R  is  the  target 
range,  hk  is  the  target  height,  and  Qt  is  an  unknown  complex  amplitude  due  to  target 
characteristics.  The  complex  amplitudes  of  the  direct  and  reflected  rays  are  simply  related 
by  a  complex  multipath  reflection  coefficient  A  at  point  0  as  shown  in  Fig.l.  Coefficient  A 
is  determined  from  the  reflection  and  specular  scattering  coefficients  and  divergence  factor. 
Polarization  enters  the  model  via  the  reflection  coefficient.  Standard  formulas  [29]  relate 
the  flat-earth  quantities  to  the  spherical-earth  quantities. 

The  geometry  for  the  radar  model  is  illustrated  in  Fig.2.  We  consider  four  separate 
propagation  paths  for  each  element  in  the  array:  direct  and  reflected  signals  going  from  the 
transmitter  to  the  target  and  direct  and  reflected  signals  returning  from  the  target  to  the 
receiving  array.  We  use  small  angle  approximations  [30]  to  derive  the  signal  received  at  the 
k‘*»  element  of  the  array: 

V]t  =  Q2  exp{-j2(R)  exp{-j({zl  +  hl)/2R}  exp{-j({zl  +  hl)/2R}  » 
[exp{jihTtZk/R)  +  A2  exp{-j(hi,Zk/R)][exp{jiheZe/R)  +  Ai  exp{-j^heZe/R)] 


(2.2) 
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Figure  I:  Geometry  for  the  beacon  model 


Transm  titer 


Figure  2:  Geometry  for  the  radar  model 


6 


if 


where  Ze  and  he  are  the  transmitter  and  target  heights  with  respect  to  the  tangent  plane  at 
the  reflection  point  Oi  for  the  transmitted  signal;  Q2  is  an  unknown  complex  amplitude  due 
to  target  characteristics;  Zk  and  kk  are  the  height  of  the  antenna  element  and  that  of 
the  target  measured  with  respect  to  the  tangent  plane  at  the  reflection  point,  O2,  for  the 
signal  returning  from  the  target;  and  At,  are  the  complex  multipath  reflection 
coefficients  at  point  Oi  and  O2  respectively.  All  the  quantities  are  corrected  for  the  earth 
curvature  and  these  corrections  are  different  for  each  antenna  element  height.  Therefore, 
we  represent  the  target  height  by  Ak- 

To  simplify  the  notation  let: 


Vk  =  b  fk{h) 

(2.3) 

where 

b=  Qt  exp(-j^R}  (beacon  model) 

(2.3a) 

or 

b=  Q2  exp{-j2^R)  exp{-j^{zl  +  hl)/2R}  » 
[exp(jiheZe/R)  +  At  exp{-j^heZe/R)] 

(radar 

model) 

(2.3b) 

and 

Jk(h)  =  exp{-j^{zi  +  hl)/8R}  [exp{j(hkZk/R)  +  A2  exp{-j^hkZk/R)] 
(same  for  both  beacon  and  radar  model) 

(2.3c) 

Assume  that  the  observation  vector  or  snapshot  s  coming  from  the  output  of  an 
array  having  K  sensors  is  given  by 


8  =  b  f(A)  +  71  (2.4) 

with  (s,  f,  »; )  6  ,  6  6  €‘*1  and  where  b  can  be  deterministic  (non-fluctuating  case)  or 

random  (fluctuating  case)  and  the  noise  vector  rj  is  assumed  to  be  stationary,  additive, 
spatially  white  and  independent  of  the  target  signals. 

3  The  RML  technique 


The  ML  estimate  of  an  unknown  parameter  vector  0  is  that  value  of  0  which 
maximizes  the  conditional  density  p(5|  (f)  of  the  observations,  or  likelihood  function.  The 
RML  technique  uses  ML  estimation  with  a  model  that  accounts  for  the  physics  of  specular 
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multipath.  Depending  on  the  degree  of  coherence  between  snapshots,  different  estimators 
are  derived.  We  consider  three  cases: 

1)  N  coherent  snapshots; 

2)  M  non— coherent  snapshots; 

3)  M  non-coherent  trains  of  Na  coherent  snapshots. 

3.1  N  coherent  snapshots: 

We  achieve  signal  coherence  by  using  a  coherent  radar  of  good  quality  in  order  to 
describe  the  relative  phase  shift  between  the  signal  component  of  the  snapshots  in  terms  of 
the  doppler  shift  caused  by  target  motion  between  signal  snapshots.  Pulse-Doppler  radars 
may  use  pulse  repetition  frequencies  as  high  as  25  kHz  in  some  iu^tances  although 
frequencies  of  10  kHz  or  less  are  more  usual.  The  degree  of  signal  coherence  in  these  radars 
is  necessarily  high  to  achieve  strong  clutter  rejection.  We  reject  clutter  by  using  Doppler 
filtering.  The  target  speed  is  obtained  by  identifying  the  particular  Doppler  filter  output 
which  gives  rise  to  a  detection. 

The  observation  model  is  given  by 

Sn  —  infn(h)  +  (3-ly 

where 


Si,:  represents  the  n‘^  snapshot; 

6n  =  6  exp(;¥>n)  with  yjn  the  known  Doppler  phase  shift  (determined  by  Doppler 
filtering). 

The  ML  estimate  of  0  is  that  value  of  0  which  maximizes  the  conditional  joint 
density  (i.e.  likelihood  function)  or  a  monotonic  function  of  it  (i.e.  log-likelihood).  The 
conditional  joint  density  is  given  by 


N 

P(8i,82,...,Sk  1  ^  CexpHSn  -  Wn)"  (Sn  “  Wn)}/a2  (3.2) 

n=l  ■ 


where  C  does  not  depend  upon  parameter  vector  0  and  N  is  the  total  number  of  snapshots. 
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The  log-likelihood  function  for  this  case,  after  dropping  the  constants  independent 

of  0,  is 

N 

=  (3-3) 

n=l 

with  0  —  [A,  6]. 

We  may  assume  that  the  model  fn(h)  and  the  noise  power  do  not  change  appreciably 
over  a  period  of  several  milliseconds  so  that  ta{h)  can  be  replaced  by  f(h)  and  by  cr^. 
The  log-dikelihood  function  becomes 

N 

(3-4) 

n=l 

with 

yn  =  Bn  exp{-jipn)  (3.4a) 

Here,  we  have  replaced  by  6  exp{j(pn). 

Now  let 

R 

8$  =  E  Vn  (3.5) 

n»l 

where  Ss  is  the  coherent  integration  of  the  signal  component  in  the  N  snapshots  Sn- 

Vector  Ss  can  always  be  formed  if  we  know  tpn.  It  is  the  a  priori  knowledge  of  the  (pn 
which  defines  the  snapshots  as  coherent.  The  equation  (3.4)  can  now  be  written  as 

m  =  -  +  N  IIHI“  Pl“]  (3.6) 

A  necessary  condition  for  L{ff)  to  be  maximized  is  obtained  by  setting  the  partial 
derivative  of  (3.6)  with  respect  to  b  [31]  to  zero.  This  condition  yields  an  estimate  of  6  as 


Substituting  (3.7)  in  (3.6),  we  obtain  the  ML  estimator  of  the  height  of  a  low-altitude 
target  as  the  value  of  h  that  maximizes  the  following  function: 


(^{h)  = 


M 


N(^S,||8n||2)||f(;i)||2 

where  the  superscript  ^  means  coherent. 


(3.8) 


We  determine  the  ML  estimate  of  h  by  adjusting  h  in  C^{h)  until  we  find  the  largest 
peak.  The  summation  of  the  squared  magnitudes  of  the  Sn  snapshots  in  the  denominator  of 
C^{h)  is  simply  a  scaling  factor  which  can  be  replaced  by  ||Ss||2/N.  In  this  case,  the 
estimator  has  exactly  the  same  form  as  obtained  for  a  single  snapshot  except  that  s  is 
replaced  by,  Ss,  the  coherent  integration  in  (3.5).  Using  the  Cauchy— Schwartz  inequality, 
we  show  that 


N  E  IlSnp  =  N  S  livnip  >  II J  YnlP  =  ||Ss||2  (3.9) 

n^i  n<t  nsi 

and  we  always  have  {  0<  C^{h)  <  1). 

In  Fig.  3  is  an  example  of  C^{h)  obtained  by  simulation  using  the  following 
parameters;  a  target  range  at  5  km,  a  target  height  of  40  m,  a  vertical  array  of  8  sensors 
having  aim  aperture  centered  at  a  height  of  9.5  m  about  water,  a  smooth  sea,  and  a 
transmitted  frequency  of  10  GHz.  We  observe  that  C^{h)  has  multiple  peaks  and  presents 
height  ambiguities.  This  arises  because  the  model  is  a  mathematical  description  of  the 
interference  pattern  resulting  from  the  summation  of  a  direct  signal  and  a  reflected  or 
image  signal.  Since  the  separation  of  target  and  its  image  is  many  wavelengths,  the  two  act 
as  a  two-element  antenna  array  and  the  .height  ambiguities  correspond  to  the  grating  lobes 
of  this  antenna.  Since  this  interference  pattern  is  being  detected  by  an  array  and  the  period 
for  the  grating  lobe  repetition  is  slightly  difierent  for  each  element  in  the  array,  the  true 
peak  is  in  theory  the  largest  but  the  neighbouring  peaks  are  so  close  and  the  decay  in 
amplitude  so  slow  that  the  amplitude  is  not  a  reliable  indicator  when  operating  at  a  single 
frequency. 
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3.2  M  non-coheient  snapshots: 

Non— coherence  between  data  snapshots  can  have  a  variety  of  causes.  One  can 
deliberately  create  non-coherence  by  the  use  of  frequency  agility;  an  important  benefit  is 
the  resolution  of  height  ambiguities  described  in  the  previous  section.  Another  cause  of 
non-coherence  is  the  large  time  separation  between  the  sampling  of  groups  of  snapshots 
sometimes  encountered  with  experimental  systems.  Still  other  causes  are  the  lack  of 
coherence  between  reference  oscillators  or  the  lack  of  knowledge  of  the  target  velocity. 

We  assume  that  the  observation  vectors  are  statistically  independent  from  snapshot 
to  snapshot.  The  noise  levels  may  differ  for  each  frequency  as  a  result  of  the  frequency 
dependence  of  the  receiver  noise  figure.  The  observation  model  is  given  by 

8in  =  6aifin(h)  +  Tjn  (3.10) 

The  log-likelihood  function  for  this  case,  after  dropping  the  constants  independent 

of  0,  is 


H 

m=l 


(3.11) 


where  0=  [h,  bi,  b2,  ...  ,6m]  and  M  is  the  number  of  snapshots  taken  at  different  frequencies. 

A  necessary  condition  for  L(0)  to  be  maximized  is  satisfied  by  setting  the  partial 
derivative  of  (3.11)  vdth  respect  to  each  bn  [31]  to  zero.  This  condition  yields  an  estimate  of 
6iii  as 


(fS  »■) 

w 


m=l,2,...,M 


(3.12) 


Substituting  (3.12)  in  (3.11),  we  obtain  the  ML  estimator  of  the  target  height  as  the  value 
of  h  that  maximizes  the  following  function: 


(h) 


Height  <m) 


Figure  3  An  example  of  C^{h) 


Height  (m) 


Figure  4:  An  example  of  0^{h) 
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\KfMP 


(3.13) 


where  the  superscript  °  means  non  coherent  snapshots.  By  using  the  Cauchy— Schwartz 
inequality,  we  can  show  that  0  <  C^{h)  <  1 . 


In  Fig.  4  is  an  example  of  C^(h)  obtained  by  simulation  using  the  same  parameters 
as  in  Fig.3  but  now  with  two  transmitted  radar  frequencies  of  9  and  10  GHz.  Using  more 
than  two  frequencies  accentuates  the  value  of  the  peak  corresponding  to  the  true  target 
height  relative  to  the  ambiguous  peaks  [18].  The  radar  transmits  a  pulse  at  one  frequency, 
changes  frequency  and  transmits  another  pulse. 


3.3  Mixed  coherent  and  non-<x>herent  snapshots: 


The  problem  of  mixed  coherent  and  non-coherent  snapshots  occurs  frequently  in 
radar  as  a  result  of  the  use  of  burst-to-burst  frequency  agility.  Here  the  radar  transmits 
bursts  of  coherent  pulses  at  one  frequency,  changes  frequency  and  transmits  another  burst 
of  pulses.  Noise  levels  may  differ  on  each  pulse  burst  as  a  result  of  a  frequency  dependence 
of  the  receiver  noise  figure. 

The  observation  model  is  given  by 


Snm  —  6nmfn«(h)  +  tjnm  (^•14) 

where  m  indicates  the  m^^^  frequency  and  n  the  n‘h  snapshot  within  a  given  coherent  burst 
and  6nn  =  bn  exp{j<pnm)-  Phase  ^nn  represents  the  doppler  shift  as  a  result  of  target  motion 
and  is  assumed  known  to  within  a  constant  Oigi.  Note  that  the  unknown  oig,  is  lumped  into 
the  complex  constant  ba-  Therefore,  there  is  coherence  within  a  burst.  The  an  are  unknown 
however,  so  that  there  is  no  coherence  between  bursts  at  different  frequencies  (different 
values  of  m).  Because  we  assume  coherence  within  the  mti»  burst,  we  can  write 
fna(A)  =  fii(^)  and  consider  a  variable  level  of  receiver  noise  represented  by  for  each 
burst.  If  ^n  is  not  known,  the  model  reduces  to  the  completely  non-coherent  case  of  the 
previous  section. 
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The  conditional  joint  density  to  be  maximized  is  p(8na|0)  with  ns=l,2,...,Na  and 
m=l,2,...,M  where  Nn  indicates  the  number  of  snapshots  in  the  m^^  burst  and  M  is  the 
number  of  different  frequencies.  The  log-likelihood  function  is  then 


K 

-  X  4  L?1  WA  -  ««» + N  iiy  >  iiy 

dis  defined  as  in  the  previous  section. 

=  noise  power  for  the  m‘^i  frequency. 

Nq 

Ssffl  =  S  Vnm  =  coherent  integration  of  the  snapshots  at  the  m<^^  frequency. 

ii>  1 


(3.15) 


Following  the  same  procedures  as  in  the  previous  section  to  eliminate  of  (3.15), 
we  obtain  the  ML  estimate  of  the  height  as  the  value  of  h  that  maximizes  the  following 
function: 


M 

=  K ‘ —  y 


II  ss.y»' 

N.»jiiy» 

/ 


(3.16) 


The  superscript means  mixed  case  of  coherent  and  non-coherent  snapshots. 


4  Formulation  of  the  Cramm— Rao  bound 


The  Cramer— Rao  bound  is  a  lower  bound  to  the  variance  of  an  imbiased  estimate. 
Most  estimators  have  biases  if  the  number  of  samples  is  small  and/or  the  signal-to— noise 
ratio  small.  The  Cramer— Rao  bound  is  nevertheless  considered  an  important  indicator  of 
performance  of  ML  estimators  because  (i)  ML  estimators  are  asymptotically  unbiased  as 
the  number  of  samples  increases  and  (ii)  the  variance  of  ML  estimators  approaches  the  CR 
bound  asymptotically  as  the  number  of  samples  increases  [32]. 

The  Cramer— Rao  bound  is  formulated  in  terms  of  the  Fisher  information  matrix  J 
[32],  the  dimension  of  which  depends  on  the  number  of  unknown  parameters  in  the 
likelihood  function.  For  example,  when  the  data  vectors  are  coherently  integrated  into  a 
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single  data  vector,  there  are  three  unknown  target  parameters,  height  h,  signal  amplitude  r, 
and  signal  phase  “ip.  J  is  then  a  3x3  matrix.  In  comparison,  the  standard  angle-of— arrival 
model  has  two  unknown  amplitudes,  two  unknown  phases,  and  two  unknown  angles;  J  is 
then  a  6*6  matrix.  The  (i,j)‘^  element  of  J  is  given  by 


Jii  =  -E 


\d%log  p{s\ff)) 


dOi  dOj 


(4.1) 


and  if  JiJ  is  the  (i,j)‘*‘  element  of  J*^  then  the  Cramer— Rao  bound  on  the  variance  of  the 
estimate  of  d\  is; 


Var(tfi-^i)  >  Jii  (4.2) 

Consider  the  general  observation  model  for  the  non-coherent  snapshots  (eq.3.10). 
The  log-likelihood  function  is  then 


M 

bn  =  r*  expijipn)  and  where  0=[h,  n,  ipu  rj,  V’2, ...  ,  ru,  ipn] 

Taking  partial  derivatives  and  the  expectation  as  indicated  in  (4.1)  leads  to  the 
following  formula  [33]  for  the  Fisher  information  matrix: 


with 


M 

J  =  2y  i-Re(HKHj  glRL*L 
n=l  “ 


'^m(^  ^ni(^ 

dOi  802 


dOi 


6 


(4.4a) 


(4.4b) 


and  where  L  is  the  number  of  real  unknown  parameters  (2M+1)  and  K  indicates  the 
number  of  sensors. 
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There  is  a  total  of  L  =  (2M+1)  unknowns:  h,  {ri,r2,...,ni},  The 

matrix  has  dimension  K*(2M+1)  and  there  are  only  three  columns  with  non-^ero 
values  because 


^/kn^n  ^/kn^n 

- - - =  0  if  m  ^  n 

dlflm  drm 


Now,  partition  J  as 


where  gut  €  IR‘*S  gam  €  Ri*2,  Gan  6  RJ*^. 
and  with 


9lm  = 


t*2 

'm 


dl 


m 


dh 


g2ia  — 


rafH  1 

■ 

- Im 

— 4 

dh  J 

rafH  1 

—  Re 

— 4 

dh 

(4.4c) 


(4.5) 


(4.6a) 


(4.5b) 
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Using  the  matrix  inversion  lemma  [34,  p.41],  the  Cramer— Rao  bound  on  the 
variance  of  the  target  height  estimate  is 


with 


(4.6a) 


(4.6b) 


We  illustrate  the  Cramer-Rao  bound  for  data  vectors  that  are  coherently  integrated 
into  a  single  data  vector  by  setting  M  =  1  in  (4.6a).  Fig.  5  shows  the  Cramer— Rao  bound 
multiplied  by  SNR  and  plotted  versus  target  height  for  a  target  at  5  km,  a  radar 
wavelength  of  3  cm  and  for  various  values  of  the  amplitude  of  the  reflection  coefficient,  a, 
and  assuming  horizontal  polarization.  The  horizontal  line  represents  the  case  where  a  =  6. 
This  case  corresponds  to  the  best  that  can  be  done  with  a  model  using  a  single  plane  wave 
(monopulse,  Fourier  beamforming).  The  limit  obtained  when  k=  0,  can  be  shown  to  be 
equivalent  to  the  Cramer-Rao  limit  of  Rife  and  Boorstyn  [35,  eq.l7]  for  the  estimation  of 
the  frequency  of  a  sinewave  of  unknown  amplitude  and  phase,  provided  we  interpret  sin  0 
as  the  spatial  frequency  and  Sxd/A  as  the  interval  between  samples.  The  results  of  figure  5 
indicate  that,  over  a  very  wide  range  of  target  heights  of  interest,  modelling  the  multipath 
by  using  two  plane  waves  is  a  very  significant  aid  in  measuring  target  height  if  the 
appropriate  models  are  used.  Indeed,  the  curves  for  a=.5  and  a=.9  cross  the  a=0  curve  only 
at  target  heights  of  less  than  5  m.  At  h=5  m,  the  angular  separation  between  the  image 
and  the  direct  path  signal  is  1/25  of  a  beamwidth  for  an  array  21  wavelengths  in  extent  as 
used  in  the  example  of  Fig.5  (8  elements  with  3A  spacing). 
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Fig.5  shows  the  effect  of  sea  state  which  is  taken  into  account  by  reducing  the 
reflection  coefGdent.  The  best  performance,  as  indicated  by  the  lowest  Cramer— Rao  bound, 
occurs  for  the  lower  sea  states  (smooth  sea  a=.9).  The  ratio  of  the  CR  bound  for  a=0  to 
that  for  a=.9  approaches  30  dB.  We  expect  a  similar  ratio  between  the  error  variances  for 
high  SNR.  As  the  smooth  reflecting  surface  becomes  rough,  the  strength  of  the  specular 
multipath,  a,  decreases.  We  conjecture  that  for  suRidently  rough  seas,  modelling  the 
multipath  may  create  more  errors  than  using  a  model  with  only  the  direct  signal.  The 
effect  of  modd  errors  is  to  produce  a  bias  on  the  estimated  target  hdght. 

5.0  Results  of  Monte-Carlo  simulations 

We  have  used  the  root— mean-square  errors,  RMSE,  as  the  measure  of  performance. 
The  results  shown  in  Fig.6  use  a  linear  array  of  K=8  equally  spaced  omnidirectional 
dements,  with  an  antenna  aperture  of  one— metre  length.  The  height  of  the  target  was 
arbitrarily  fixed  at  10  metres  and  two  X— band  frequendes  were  used  (9  and  10  GHz) 
operating  over  a  smooth  sea  corresponding  to  Sea  State  2  or  lower.  We  carried  out 
Monte-Carlo  simulations  to  determine  the  threshold  effects  on  height  estimation.  We  have 
also  plotted  the  Cramer-Rao  bound.  Our  estimators  lie  very  close  to  the  Cramer-Rao 
limit  for  large  SNR  but  large  deviations  occur  as  the  SNR  decreases  bdow  a  threshold 
value. 


The  term  SNR  ,  is  a  per— element  averaged  quantity  summed  over  the  number  of 
frequencies: 


SNRdB  =  10  log 


2K  al 


(5.1) 


where  (M=2)  is  the  number  of  frequencies,  ht=10  m  is  the  target  hdght  and  is  the  noise 
power. 


We  have  plotted  on  Fig.6,  the  Cramer-Rao  bound  and  the  estimator  0{h).  Our 
estimator  lies  very  close  to  the  Cramer-Rao  limit  for  large  SNR  but  large  deviations  occur 
as  the  SNR  decreases  below  a  threshold  value.  To  decrease  this  estimation  threshold,  we 
increase  the  number  of  frequencies  and  the  frequency  bandwidth. 
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Figure  5:  Comparison  of  CR  bounds  with  various  values 
of  reflection  coefficient 
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Figure  6:  Threshold  effects 
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6.  EXPERIMENTAL  RESULTS. 

The  performance  of  the  RML  technique  has  been  evaluated  with  two  experimental 
systems:  the  ELAT  and  the  MARS  systems.  The  ELAT  results  are  presented  in 
Figs.  7  to  11.  Each  figure  has  three  parts:  1)  estimated  SNR,  2)  FFT  beamforming  results, 
3)  the  RML  results. 

We  estimate  the  noise  power  by  using  a  projection  technique  to  remove  the  specular 
component  as 


9 


(IS(At)s.) 

nm  —  8|b - ; - fffl(ht) 

WUihW 

where  ht  is  the  ML  estimate  of  the  target  height  obtained  from  our  analysis  of  C^{h).  The 

quantity  ||na,|p  is  an  estimate  of  the  combined  power  of  the  receiver  noise  and  the  diffuse 
multipath  for  the  mt^i  frequency. 

The  estimated  SNR  is  then  given  by  the  following  equation: 


SNR(dB)  =  10  log 


llnalP  1 


(6.2) 


The  estimate  given  by  (6.2)  is  in  reality  the  ML  estimation  of  SNR.  By  the  invariance 

property  of  ML  estimators  [36,  p.92],  we  have  that  g(^L)  is  the  ML  estimate  of  g(0).  If  we 
apply  this  property  to  our  case  the  ML  estimate  of  the  mean  6nifin(At)  is  given  by 


Ab  —  8(^l)  ~  6mfm(ht)  (6‘3) 

with  hm  given  by  the  relation  (3.12). 

It  is  well  known  that  the  ML  estimate  of  the  mean  of  a  Gaussian  population  is  equal 
to  the  sample  mean  and  the  ML  estimate  of  the  variance  of  a  Gaussian  population  is 
simply  equal  to  the  sample  variance  [36,  p.91].  The  noise  variance  is  then  given  by 
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a2=||8.-;in.||2  (6.4) 

When  we  compare  (6.4)  with  the  amplitude  squared  of  (6.1),  we  conclude  that  the 
projection  method  yields  the  ML  estimate  of  the  noise  power 

The  middle  part  on  each  figure  presents  results  obtained  with  the  classical  Fourier 
beamforming.  We  choose  to  compare  our  results  with  those  of  Fourier  beamforming 
because  Fourier  beamforming  produces  angular  tracking  errors  comparable  to  monopulse, 
the  most  vndely  used  tracking  radar.  It  is  anticipated  that  superresolution  algorithms 
[6—7],  when  applicable  to  low— angle  tracking,  produces  angular  tracking  errors  between 
those  obtained  with  Fourier  beamforming  and  the  RML  technique.  This  topic  will  be 
covered  in  a  subsequent  paper. 

The  results  for  conventional  beamforming  are  obtained  by  forming  beams  directly  to 
obtain  an  estimate  of  the  target  height.  The  set  of  the  eight  sensor  outputs  collected  at  a 
particular  instant  of  time  are  augmented  with  zeros  to  form  a  complex  vector  of  1024 
points.  We  use  the  Fast  Fourier  Transform  (FFT)  in  the  angular  domain  to  obtain  an 
estimate  of  the  elevation  angle,  which  is  converted  tc  a  height  estimate.  This  is  repeated 
for  the  second  operating  frequency  and  the  final  estimate  is  an  average  of  the  two 
estimates.  This  feature  significantly  improves  the  estimation  performance  compared  with 
single-frequency  conventional  beamforming;  the  results  of  which  are  not  shown  here. 

The  bottom  part  of  each  figure  presents  results  obtained  using  the  RML  technique 
with  the  estimator  C'*(h). 

6.1  Experimental  results  using  the  ELAT  system: 

This  section  presents  examples  of  experimental  target  tracking  using  the  ELAT 
radar.  ELAT  is  being  developed  at  DREO  as  part  of  a  research  program  into  high  angular 
resolution  techniques  for  radar.  The  ELAT  radar  is  designed  with  a  sampled  aperture 
antenna  composed  of  8  receiving  elements  spaced  at  a  distance  of  14  cm.  Each  element  is 
coimected  to  an  individual  receiver  for  conversion  of  8.65  or  9.6  GHz  microwave  frequency 
information  to  the  video  frequency  band.  Each  receiver  includes  quadrature  (Q)  and 
in— phase  (I)  detectors  which  are  followed  by  analog— to-digital  conversion  in  each  channel. 
Using  this  receiving  arrangement,  discrete  samples  of  the  electromagnetic  field  across  the 
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Figure  7:  Tracking  an  ascending  target  over  smooth  sea  (SSI) 
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Figure  8:  Tracking  an  ascending  target  over  smooth  sea  (SSl) 
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Figure  9:  Tracking  an  ascending  target  over  smooth  sea  (SSI) 
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vertical  antenna  aperture  can  be  converted  to  the  video  band  to  form  the  observation 
vector  or  snapshot  s,  used  in  the  previous  sections. 

To  test  C“(h)  we  use  five  files  collected  on  the  Ottawa  River.  We  have  installed  a 
beacon  in  a  low-flying  helicopter.  The  two  frequencies  used  are  8.65  and  9.6  GHz.  For  the 
first  three  examples  the  helicopter  moves  vertically  from  approximately  10  m  to  350  m. 
This  profile  is  repeated  for  three  different  ranges  (4.7  km,  6.7  km,  9.7  km).  The  other  two 
files  represent  a  target  manoeuvring  in  height  over  a  15  km  path  with  the  helicopter  flying 
toward  the  vertical  array.  The  solid  line  on  each  figure  represents  the  theoretical  pattern 
that  the  pilot  tried  to  follow. 

From  Figs. 7— 9,  we  can  make  the  following  observations: 

1.  The  variance  on  the  height  estimation  increases  with  range  (Fig.9  compared  to 

Fig.7); 

2.  The  RML  method  using  C"(/i)  gives  superior  results  when  compared  to  Fourier 
beamforming  specially  at  low  altitude,  and 

3.  Fourier  beamforming  gives  useful  estimates  only  when  the  angular  separation 
between  the  target  and  its  image  is  one  beamwidth  or  larger.  In  this  example,  one 
beamwidth  is  approximately  equal  to  0.034  (=A/D,  where  D  is  the  array  aperture) 
radians,  corresponding  to  a  target  height  of  80  m  at  4.7  km  or  160  m  at  9.7  km. 

The  next  two  figures.  Figs. 10— 11,  present  the  tracking  over  a  15  km  path.  The 
results  indicate  that  our  technique  gives  an  useful  estimate  of  the  target  height  while  the 
classical  Fourier  beamforming  produced  poor  target  height  estimates.  The  results  obtained 
with  the  RML  technique  are  relatively  good  considering  that  the  data  have  been  collected 
on  the  Ottawa  River  over  a  propagation  path  that  was  probably  non— homogeneous  due  to 
the  proximity  of  land. 

6.2  Experimental  results  using  the  MARS  system 

The  performance  of  the  RML  algorithm  increases  when  using  a  larger  bandwidth 
and  multiple  frequencies.  This  is  in  agreement  with  the  experiment  results  derived  from 
data  collected  by  McMaster  University  at  a  site  on  Lake  Huron  using  an  array  system 
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Figure  11:  Tracking  a  target  maneuvering  over  smooth  sea  (SSI) 
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called  MARS.  The  MARS  system  comprises  a  multiple  frequency  beacon  system,  a 
sampled— aperture  antenna  and  a  data  acquisition  system.  Two  frequencies  are 

*  simultaneously  used  for  transmission:  a  fixed  frequency  of  10.2  GHz  and  an  agile  frequency 
ranging  Rom  8.02  to  12.4  GHz  in  30— MHz  steps.  In  the  experiments,  only  18  frequencies 

1.  were  actually  used.  The  sampled  aperture  contains  a  linear  array  of  32  elements  uniformly 

spaced  over  a  1.82  m  aperture.  Each  antenna  element  has  its  own  receiver  and  the  data  are 
digitized  with  12— bit  precision. 

For  the  Lake  Huron  experiments  [12—13],  the  transmitter  was  sited  at  a  distance  of 
4610  m  from  the  receiving  antenna.  Fig.l2  is  a  picture  of  the  measurement  setup.  The 
height  of  the  transmitter  can  vary  from  2.7  m  to  18.6  m  above  the  water  level.  The  height 
of  the  waves  was  visually  observed  to  be  about  1.5  m  to  2  m  (valley  to  peak)  corresponding 
to  a  sea  state  of  3. 

In  Fig.l3  we  present  the  results  obtained  with  a  target  varying  its  height  Rom  3.5  m 
(.1  beamwidth)  to  15.5  m  (.4  beamwidth).  The  target  range  is  4610.4  m.  The  theoretical 
values  are  indicated  by  the  solid  line  in  the  figure  while  the  estimated  target  heights  are 
indicated  by  a  small  triangle.  We  observe  that  the  variance  increases  when  the  angular 
separation  between  the  target  and  its  image  decreases.  Biases  are  observed  in  all  results. 
The  results  are  obtained  over  a  moderately  rough  sea  surface  (SS3)  characterized  by  visual 

*  information  and  bias  is  created  if  the  sea  state  is  over  or  under  estimated.  Finally  tilt  and 
vibrations  caiised  by  strong  wind  (18  km/hr)  contribute  to  both  bias  and  the  variance  on 
the  estimation.  In  Fig.  14,  both  the  variance  and  the  bias  decrease  when  we  increase  the 
bandwidth  and  the  number  of  Requencies.  The  SNR  estimated  for  the  examples  presented 
in  Figs.l3— 14  is  approximately  30  dB. 

7.  CONCLUSION 

We  have  presented  the  theory  for  a  new  form  of  maximum  likelihood  estimation  for 
low— angle  tracking.  This  method  is  based  on  a  "refined"  propagation  model  employing  a 
priori  information;  we  have,  therefore,  coined  the  name  Refined  Maximum  Likelihood 
(RML).  The  results  of  simulations  and  experiments  with  beacon  signals  indicate  substantial 

'  improvements  compared  with  previous  techniques.  However,  the  RML  model  is  imperfect; 

it  does  not  include  diffuse  multipath,  sea  swell  or  ducting.  The  results  of  limited 

*  experiments  with  beacons  systems  in  estimated  sea  state  3  conditions  are  encouraging  and 
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indicate  some  robustness  in  the  presence  of  diffuse  multipath.  More  extensive 
experimentation  with  radar  signals  is  required  to  corroborate  these  results  and  to 
determine  the  applicability  of  the  RML  and  the  limitations  of  this  applicability  in  realistic 
environments.  Finally  we  want  to  emphasize  the  benefits  of  wide  band  frequency  agility  for 
effective  resolution  of  ambiguities. 
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Figure  12:  Geometry  of  the  Lake  Huron  experimental  setup 
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Figure  13:  Example  of  tracking  over  relatively  rough  sea  (SS3) 
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Figure  14:  Effects  of  frequency  agility  on  tracking  (SS3) 
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